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Abstract 



We propose to study the non-equilibrium features of heavy-ion collisions by 
following the evolution of an initial state with a large number of quanta with 
a distribution around a momentum |fco| corresponding to a thin spherical shell 
■ in momentum space, a 'tsunami'. An O(N) (<I> 2 ) 2 model field theory in the 

large N limit is used as a framework to study the non-perturbative aspects of 
Qh! the non-equilibrium dynamics including a resummation of the effects of the 

Q_i! medium (the initial particle distribution). In a theory where the symmetry is 

spontaneously broken in the absence of the medium, when the initial number 
of particles per correlation volume is chosen to be larger than a critical value 
the medium effects can restore the symmetry of the initial state. We show 
that if one begins with such a symmetry-restored, non-thermal, initial state, 
non-perturbative effects automatically induce spinodal instabilities leading to 
a dynamical breaking of the symmetry. As a result there is explosive particle 
production and a redistribution of the particles towards low momentum due to 
the nonlinearity of the dynamics. The asymptotic behavior displays the onset 
of Bose condensation of pions and the equation of state at long times is that of 
an ultrarelativistic gas although the momentum distribution is non-thermal. 

11.10.-z,11.30.Qc,11.15.Tk 
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I. INTRODUCTION 



The Relativistic Heavy Ion Collider (RHIC) at Brookhaven and the Large Hadron Collider 
(LHC) at CERN will provide an unprecedented range of energies and luminosities that will 
hopefully probe the Quark-Gluon Plasma and Chiral Phase transitions. The basic picture of 
the ion-ion collisions in the energy ranges probed by these accelerators as seen in the center- 
of-mass frame (cm.), is that of two highly Lorentz-contracted 'pancakes' colliding and leaving 
a 'hot' region at mid-rapidity with a high multiplicity of secondaries [ffj. At RHIC for Au+Au 
central collisions with typical luminosity of 10 26 /cm 2 . s and cm. energy pa 200GeV/n — n, a 
multiplicity of 500-1500 particles per unit rapidity in the central rapidity region is expected 
Q-fp. At LHC for head on Pb + Pb collisions with luminosity 10 27 /cm 2 . s at cm. energy 
pa 5 TeV/n — n, the multiplicity of charged secondaries will be in the range 2000 — 8000 per 
unit rapidity in the central region At RHIC and LHC typical estimates |l]-§| of energy 
densities and temperatures near the central rapidity region are e ~ 1 — 10 GeV/fm 3 , T ~ 
300- 900MeV. 

Since the lattice estimates |||| of the transition temperatures in QCD, both for the QGP 
and Chiral phase transitions are T c pa 160 — 200 MeV, after the collision the central region 
will be at a temperature T > T c . In the usual dynamical scenario that one |l| envisages, the 
initial state cools off via hydrodynamic expansion through the phase transition down to a 
freeze-out temperature, estimated to be Tp pa 100 MeV 0, at which the mean free-path of 
the hadrons is comparable to the size of the expanding system. 

The initial state after the collision is strongly out of equilibrium and there are very few 
quantitative models to study its subsequent evolution. There are perturbative and non- 
perturbative phenomena that contribute to the processes of thermalization and hadroniza- 
tion. The perturbative (hard and semihard) aspects are studied via parton cascade models 
which assume that at large energies the nuclei can be resolved into their partonic constituents 
and the dynamical evolution can therefore be tracked by following the parton distribution 
functions through the perturbative parton-parton interactions []8}-|i"2"[ . Parton cascade models 
(including screening corrections to the QCD parton-parton cross sections) predict that ther- 
malization occurs on time scales ~ 0.5fm/c [fL3|l . After thermalization, and provided that the 
mean-free path is much shorter than the typical interparticle separation, further evolution 
of the plasma can be described with boost-invariant relativistic hydrodynamics The 



details of the dynamical evolution between the parton cascade through hadronization, and 
eventual description via hydrodynamics is far from clear but will require a non-perturbative 
treatment. The non-perturbative aspects of particle production and hadronization typically 
envisage a flux-tube of strong color-electric fields, in which the field energy leads to produc- 
tion of qq pairs [|IJ],|llJ. Recently the phenomenon of pair production from strong electric 
fields in boost-invariant coordinates was studied via non-perturbative methods that address 
the non-equilibrium aspects and allow a comparison with hydrodynamics ]17|] . 

The dynamics near the phase transition is even less understood and involves physics be- 
yond the realm of perturbative methods. For instance, considerable interest has been sparked 
recently by the possibility that disoriented chiral condensates (DCC's) could form during the 
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evolution of the QCD plasma through the chiral phase transition [|T^]- P3] . Rajagopal and 
Wilczek |24] have argued that if the chiral phase transition occurs strongly out of equilibrium, 
spinodal instabilities [ES[] could lead to the formation and relaxation of large pion domains. 
This phenomenon could provide a striking experimental signature of the chiral phase transi- 
tion and could provide an explanation for the Centauro and anti-Centauro (JACEE) cosmic 
ray events [[2f|. An experimental program is underway at Fermilab to search for candidate 
events [p7t p8||. Most of the theoretical studies of the dynamics of the chiral phase transition 
and the possibility of formation of DCC's have focused on initial states that are in local 
thermodynamic equilibrium (LTE) ||29|| - [32 . 



We propose to study the non- equilibrium aspects of the dynamical evolution of highly 
excited initial states by relaxing the assumption of initial LTE (as would be appropriate 
for the initial conditions in a heavy-ion collision). Consider, for example, a situation where 
the relevant quantum field theory is prepared in an initial state with a particle distribution 
sharply peaked in momentum space around fco and — ko where kg is a particular momentum. 
This configuration would be envisaged to describe two 'pancakes' or 'walls' of quanta moving 
in opposite directions with momentum \k \. In the target frame this field configuration 
would be seen as a 'wall' of quanta moving towards the target and hence the name 'tsunami' 
33fl . Such an initial state is out of equilibrium and under time evolution with the proper 



interacting Hamiltonian, non-linear effects should result in a redistribution of particles, as 
well as particle production and relaxation. The evolution of this strongly out of equilibrium 
initial state would be relevant for understanding phenomena such as formation and relaxation 
of chiral condensates. Starting from such a state and following the complete evolution of the 
system thereon, is clearly a formidable problem even within the framework of an effective 
field theory such as the linear cr-model. 

In this article we consider an even more simplistic initial condition, where the occupation 
number of particles is sharply localized in a thin spherical shell in momentum space around 
a momemtum |/c |, i.e. a spherically symmetric version of the 'tsunami' configuration. The 
reason for the simplification is purely technical since spherical symmetry can be used to 
reduce the number of equations. Although this is a simplification of the idealized problem, 
it will be seen below that the features of the dynamics contain the essential ingredients to 
help us gain some understanding of more realistic situations. 

We consider a weakly coupled A$ 4 theory (A ~ 10~ 2 ) with the fields in the vector 
representation of the O(N) group. Anticipating non-perturbative physics, we study the 
dynamics consistently in the leading order in the 1/N expansion which will allow an analytic 
treatment as well as a numerical analysis of the dynamics. 

The pion wall scenario described above is realized by considering an initial state described 
by a Gaussian wave functional with a large number of particles at \ko\ and a high density is 
achieved by taking the number of particles per correlation volume to be very large. As in 
finite temperature field theory, a resummation along the lines of the Braaten and Pisarski 
|34j program must be implemented to take into account the non-perturbative aspects of the 
physics in the dense medium. As will be explicitly shown below, the large N limit in the case 
under consideration provides a resummation scheme akin to the hard thermal loop program 
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[34]. 

The dynamical evolution of this spherically symmetric "tsunami" configuration described 
above reveals many remarkable features: i) In a theory where the symmetry is spontaneously 
broken in the absence of a medium, when the initial state is the O(N) symmetric, high densty, 
"tsunami" configuration, we find that there exists a critical density of particles depending 
on the effective (HTL-resummed) coupling beyond which spinodal instabilities are induced 
leading to a dynamical symmetry breaking, ii) When these instabilities occur, there is 
profuse production of low-momentum pions (Goldstone bosons) accompanied by a dramatic 
re-arrangement of the particle distribution towards low momenta. This distribution is non- 
thermal and its asymptotic behavior signals the onset of Bose condensation of pions. iii) 
The final equation of state of the "pion gas" asymptotically at long times is ultrarelativistic 
despite the non-equilibrium distributions. 

The paper is organized as follows: In Section II we introduce the model under consid- 
eration and describe the non-perturbative framework, namely the large N approximation. 
Section III is devoted to the construction of the wave functional and a detailed description 
of the initial conditions for the problem. The dynamical aspects are covered in Section IV. 
We first outline some issues dealing with renormalization and then provide a qualitative 
understanding of the time evolution using wave functional arguments. We argue that the 
system could undergo dynamical symmetry breakdown and provide analytic estimates for 
the onset of instabilities. We present the results of our numerical calculations in Section 
V which confirm the robust features of the analytic estimates for a range of parameters. 
In Section VI we analyze the details of symmetry breaking and argue that the long time 
dynamics can be interpreted as the onset of formation of a Bose condensate even when the 
order parameter vanishes. 

Finally in Section VII we present our conclusions and future avenues of study. 



II. THE MODEL 

As mentioned in the introduction we consider a A$ 4 theory with O(N) symmetry in the 
large- N limit with the Lagrangian, 

C = - ■ $) - A ( $ . $f ( 2.i) 

where $ is an O(N) vector, $ = (cr, 7?) and 7? represents the N—l pions, 7? = (n 1 , n 2 , n 1 ^^ 1 ). 
We then shift a by its expectation value in the non-equilibrium state 



a(x, t) = VN<f>(t) + x(x, t) ; (a(x, t)) = VN<f>(t) . (2.2) 
We refer the interested reader to several articles which discuss the implementation of the 



large N limit (see for e.g. p5|-p9|, ^4] , ^5|1 ) . The 1/N series may be generated by introducing 
an auxiliary field a(x) which is an algebraic function of $ 2 (x), and then performing the 
functional integral over a(x) using the saddle point approximation in the large N limit 



4 



|35H37| . It can be shown that the leading order terms in the expansion can be easily obtained 
by the following Hartree factorisation of the quantum fields ]38|39|,|4j,45|, 



X 



6(X 2 )X 2 + constant ; x 3 -> Hx 2 )x 



(7T • T?) 2 -> 2(7f 2 )7T i - (7T 2 ) 2 + <D(1/N) 



^2 2 



^ 2 <x 2 > + (^ 2 )x 2 ; ^ 2 x - (v? 2 )x 



(2.3) 



All expectation values are to be computed in the non-equilibrium state. In the leading order 
large N limit we then obtain, 



C 



\* ■ («9 2 + Ml(t))n - \ X {d 2 + Ml(t)) X - xV'mit) + ^V> 2 



Ml(t) =m% + 



Ml{t) 



v'W),t) 

(- 2 ) 



A r 

2 
A 



m 2 R H — 
2 



<P\t) + (vr 2 ) 
30 2 (t) + (vr 2 ) 

(vr 2 )/iV. 



(2.4) 
(2.5) 
(2.6) 

(2.7) 
(2.8) 



This approximation allows us to expand about field configurations that are far from the 
perturbative vacuum. In particular it is an excellent tool for studying the behaviour of 
matter in extreme conditions such as high temperature or high density |1T7| , |35H39|j4^j45[1 . 

One way to obtain the non-equilibrium equations of motion is through the Schwinger- 
Keldysh Closed Time Path formalism. This is the usual Feynman path integral defined on a 
complex time contour which allows the computation of in-in expectation values as opposed 
to in-out S-matrix elements. For details see the references | 40| . The Lagrangian density in 
this formalism is given by 



C neq = £[$+] - £[$" 



(2.9) 



with the fields $ ± defined along the forward (+) and backward (— ) time branches. The 
non-equilibrium equations of motion are then obtained by requiring the expectation value 
of the quantum fluctuations in the non-equilibrium state to vanish i.e. from the tadpole 
equations [|3Tfl , 



(2.10) 



In the leading order approximation of the large iV limit, the Lagrangian for the x field is 
quadratic plus linear and the tadpole equation for the x leads to the equation of motion for 
the order parameter or the zero mode 



+ ^[0 2 (t) + (7r 2 )(t)]0(t)+m|0(t) 



. 



(2.11) 
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In this leading approximation the non-equilibrium action for the N — 1 pions is, 



J rf 4 x[£ w+ - C n -} = J d 3 xdt • <9 2 7f+ - Mf(t)n + ■ 7f+) - (+ — > -)} . (2.12) 

We have not written the action for the x field fluctuations because they decouple from the 
dynamics of the pions in the leading order in the large N limit [38,[39|. 

Having introduced the model and the non-perturbative approximation scheme the next 
step is to construct an appropriate non-equilibrium initial state or density matrix. 

Although one could continue the analysis of the dynamics using the Schwinger-Keldysh 
method, we will study the dynamics in the Schrodinger representation in terms of wave- 
functionals because this will display the nature of the quantum states more clearly. We find 
it convenient to work with the Fourier-transformed fields defined as, 

^*) = 4rE e "* af ^(*). ( 2 - 13 ) 

VV k 

where we have chosen to quantize in a box of finite volume V that will be taken to infinity 
at the end of our calculations. The Hamiltonian for the pions is given by 

= £ (% ■ n_* + Um ■ u) - ^ (Uvk ■ U)) . ( 2 - 14 ) 

where 



2 k -k 2 

k " \ k 



4(t) = P + Ml(t) (2.15) 

is the effective time dependent frequency and M^{t) is given by Eq. Q2.5Q . To leading order 
in the large N limit the theory becomes Gaussian and the non-linearities are encoded in a 
self-consistency condition, since the frequency Q2.15Q depends on (it 2 ) and this expectation 
value is in the time dependent state, as displayed by the set of equations (|2.5|-|2~8|). 



III. THE INITIAL STATE 



As stated in the introduction, our ultimate goal is to model and study the non-equilibrium 
aspects of the evolution of an initial, highly excited state that relaxes following high energy, 
heavy-ion collisions. An idealized description of the associated physics would be to consider 
two wave packets made up of very high energy components representing the heavy ions and 
moving with a highly relativistic momentum toward each other. The goal would be to follow 
the dynamical evolution of the wavefunctionals corresponding to this situation, thus clearly 
elucidating the non-equilibrium features involved in the phase transition processes following 
the interactions of the wave packets. This initial state could be described by a distribution 
of particles, sharply peaked around some special values ko and — ko in momentum space. 
The evolution of this state then follows from the functional Schrodinger equation. 
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Even with the simplification of a scalar field theory such a program is very ambitious 
and beyond the present numerical capabilities. One of the major difficulties is that selecting 
one particular momentum breaks rotational invariance and the evolution equations depend 
on the direction of wave vectors even in the Gaussian approximation. (This statement will 
become clear below). 

In this article however, we choose to study a much simpler description of the initial state 
which is characterized by a high density particle distribution in a thin spherical 'shell' in 
momentum space . We propose an initial particle distribution that has support concentrated 
at \ko\. This particular state does not provide the necessary geometry for a heavy ion collision, 
however it does describe a situation in which initially there is a large multiplicity of particles 
in a small momentum 'shell', there is no special beam- axis and the pions are distributed 
equally in all directions with a sharp spatial momentum. This is a rotation invariant state 
that describes a highly out of equilibrium situation and that will relax during its time 
evolution (a spherical "tsunami"). 



A. The Wave Functional: 



Since in the leading order approximation in the large N expansion the theory has become 
Gaussian (at the expense of a self-consistency condition), we choose a Gaussian ansatz for 
the wave-functional at t — 0. The reason for this choice is that upon time evolution this 
wave functional will remain Gaussian and will be identified with a squeezed state functional 
of pions. 



tf(t = 0) =%A4(0)exp 



A fc (0) 



(3.1; 



This state will then evolve according to the Hamiltonian given in Eq. ( |2.14| ) which is es- 
sentially a harmonic oscillator Hamiltonian with self-consistent, time-dependent frequencies. 
The functional Schrodinger equation is given by 



TT T 



(3.2) 



The last term in the Hamiltonian ( |2.14| ) which is independent of the fields (a time dependent 
'vacuum energy term') can be absorbed in an overall time dependent phase of the wave 
functional. Removing this term by a phase redefinition, the functional Schrodinger equation 
becomes 



ztf fo] = J2 

k 



h 2 5 2 



2 S%5ff_ % 



(3.3) 



which then leads to a set of differential equations for the covariance Af.. The time dependence 
of the normalization factors A4 is completely determined by that of the Ak as a consequence 
of unitary time evolution. The state for arbitrary time t takes then the form: 
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tf(t) = %jV fc (t)exp 



A fc (t) 



(3.4) 



The evolution equations for the covariance are obtained by taking the functional deriva- 
tives and comparing powers of on both sides. We obtain the following evolution equations 



32.38 



iA k {t) = Al(t)-c4(t), 

A4(t) = A4(0) exp [ / A Ik (t')dt' 



(3.5) 
(3.6) 



with Ak = Ankif) + iAjkit). The equal time two-point correlation function in the time 
evolved non-equilibrium state is given by 



<*\fj n -fi_ n \v> 

< > 



IWi (Vk-v-k) nA/V exp 


' AM) - - ' 
t Vq ■ V-q\ 


J[Vff ? ] TIM? exp 


" AM) -> - 

2 Vq- 7 l- q _ 





N 



2A Rk (t) 

leading to the self-consistency condition 

<- 2 >W = E 



(3.7) 



2A m (t) 



(3.1 



Formally, one can also represent these two-point equal time correlators in terms of func- 
tional integrals over the closed time path contour where the initial state is chosen to be 
the Gaussian functional described above. However the explicit and rather simple ansatz 
for the wave functional enables one to obtain the two-point functions directly in a rather 
straightforward manner. Moreover, the wave functional approach will permit a much clearer 
understanding of the physics of the problem. The Ricatti equation ( p.5| ) can be cast in a 
simpler form by writing Af. in terms of new variables <pl as 



A k {t) 



leading to the simple equation for the new variables 

r k + u 2 k (t) <i>i = o. 



(3.9) 



(3.10) 



In terms of these mode functions we find that the real and imaginary parts of the covariance 
Ak are given by 
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A Ik (t) 



d_ 

"dt 



\n\Mt)\' 



From the differential equation for the <f>k(t) given by Eq. 



(3.12) 

TUf) it is clear that the 



combination that appears in the numerator of Eq. (|3 . 1 1 ) is the Wronskian Qf. of the differ- 
ential equations and will consequently be determined from the initial conditions alone. The 
expression for the quantum fluctuations (n 2 ) = (ti 2 )/N is given by, 



d 3 k 



d 3 k |0 fc (t)p 
(2tt) 3 2VL k 



(3.13) 



The mode functions <pk have a very simple interpretation: they obey the Heisenberg 
equations of motion for the pion fields obtained from the Hamiltonian ( |2.14[) . Therefore we 
can write the Heisenberg field operators as 



jr(x, t) 



E 



4 (t>k{t) e %k - s + a\ <f>* k {t) e 



-ik-x 



(3.14) 



where a k , a\ are the time independent annihilation and creation operators with the usual 
Bose commutation relations. 



B. Initial Conditions: 



Within this Gaussian ansatz for the wave functional, the initial conditions are completely 
determined by the initial conditions on the mode functions (j>k{t). In order to physically 
motivate the initial conditions we now establish the relation between the particle number 
distribution and these mode functions. 

Since, in a time dependent situation there is an ambiguity in the definition of the particle 
number, we define the particle number with respect to the eigenstates of the instantaneous 
Hamiltonian ( |2.14|) at the initial time, i.e. 



1 5 2 



cj fc (0) V 2 5%^ 
1 



+ ^P-%-V-k 



Wfe(O) 



(3.15) 



Here, cufc(0) is the frequency (j2.15|) evaluated at t = 0, i.e. the curvature of the potential 
term in the functional Schrodinger equation ( [3.3|) at t = and provides a definition of the 
particle number (assuming that u>l(0) > 0). The expectation value of the number operator 
in the time evolved state is then 
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\A Rk (t) - uj k (0)} 2 + A 2 Th (t) 

A fc (t) 2 + ^(t) 2 

~ 4[1 + A*(t)] ' (3 - 17) 
where A k and 5 k are defined through the relations, 

A Rk (t) = u k {0) [1 + A k (t)] ; A Ik (t) = cu fc (0) <f fc (t) . (3.18) 

In terms of the mode functions (f) k and 4> k the expectation value of the number operator is 
given by 

"*(*) = I o 1 ,nx [|4W| 2 + ^ 2 (O)|0 fc (t)| 2 l - ~ . (3.19) 



4 fi fc w fc (0) 



2 



The quantity 5^(t) appears as the phase of the wave function and will be chosen to be zero 
at t — for simplicity, 

A Ik (0) = ; 4(0) = 0. (3.20) 

Assuming 5 k (0) = 0, the initial conditions on the <f) k (0) variables can be obtained at once 
from Eq. ( |3.9| ) and are found to be, 

fa(0)=in k = iu k (p)[l + A k (0)] ;<p* k (0) = l, (3.21) 

where fl k is the Wronskian W[<j) k (t)* ,4> k (t)]. Hence the wave functional of the system at 
t — can be specified completely (up to a phase) by the single function A fc (0). 

Using Eq. ( |3.17|) one can easily solve for A^(0) in terms of the initial particle spectrum 



A k = A*(0) = 2[n*(0) ± v/n fc (0) 2 + n fc (0)] . (3.22) 

Which of the two solutions will give us interesting physics is a more subtle question that we 
shall address in the next section when we discuss the dynamics of the problem. 

Before moving on to the description of the dynamics let us briefly summarize what 
we have done. We proposed a rather simple description of a large multiplicity, high energy 
particle collision process by preparing an initial state with an extremely high number density 
of particles concentrated at momenta given by \k\ = k . Consistent with the leading order in 
a large iV approximation, we chose a Gaussian ansatz for our wave functional, parametrized 
by the variables 4> k (t), <fi k (t) (or alternatively A k (t) and 5 k (t)) and the initial conditions on 
these variables are determined completely by the choice of the particle distribution function 
n k (0) at t — (Eq. ( 3.22|) ). The next step will be to obtain the renormalized equations of 
motion and then to study the dynamics analytically as well as numerically. 
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IV. THE DYNAMICS 



From the discussion in the previous sections we see that the following set of equations 
for the order parameter <f>(t) and the mode functions must be solved self-consistently 
in order to study the dynamics: 

+ + <Ob]) = > 

+ (** + ml + \[<p\t) + (tt 2 )*]) m = , (4.2) 

<P* k (0) = l ■A* k (0) = tQ k , (4.3) 
with the self-consistent condition 

{nmB = J j^r^r- (4 - 4) 

The quantities in the above equations must be renormalized. This is achieved by first 
demanding that all the equations of motion be finite and then absorbing the divergent pieces 
into a redefinition of the mass and coupling constant respectively, 

ml + ^[< vr 2 > (t) B + 2 (t)] = ml + ^[< tt 2 > (t) R + <p\t)\ = M 2 Rn (t) . (4.5) 

A detailed derivation of the renormalization prescriptions requires a WKB analysis of the 
mode functions <pk(t) that reveals their ultraviolet properties. Such an analysis has been 
performed elsewhere p8|j39|] . In summary the mass term will absorb quadratic and log- 
arithmic divergences while the coupling constant will acquire a logarithmically divergent 
renormalization |58,3~5|. In particular 



with k an arbitrary renormalization scale. 

Introducing the effective mass of the particles at the initial time as 

M 2 R = M 2 R7V (t = 0), (4.7) 

we recognize that this effective mass has contributions from the non-equilibrium particle 
distribution and is the analog of the hard-thermal loop (HTL) resummed effective mass 
in a scalar field theory. Recall, however, that the initial distribution is not thermal. In 
a scalar theory, the HTL effective mass is obtained by summing the daisy and superdaisy 
diagrams pl| which is precisely the resummation implied in the leading order in the large 



iV approximation. To see this more clearly consider the case in which the order parameter 
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vanishes, i.e. <p(t) = 0; then the effective mass at the initial time can be written as a gap 
equation 



M 2 R = m 2 B + 



m| + 



47T 5 



47T 5 



k 2 dk 

2Q k {0) 
k 2 dk 



(4- 



+ 



A 



2u k (0) An 2 J 2co k (0) 



k 2 dk 



A fc (0) 



1+A*(0) 



where we have used the relation Q k (0) = u>k(0)(l + A k ). The second term in the above 
expression is the usual contribution obtained at zero temperature (and zero density) for the 
(self-consistent) renormalized mass parameter M R i.e the one loop tadpole (with oo k (0) = 

\Jk 2 + Mr). The third term contains the non-equilibrium effects associated with the particle 
distributions and vanishes when n k — > 0. This term is finite (since n k (0) is assumed to be 
localized within a small range of momenta). For a given distribution n k (0), the solution to 
the self-consistent gap equation (|4.8|) gives the effective mass, dressed by the medium effects. 



This is indeed very similar to the finite temperature case in which the tadpole term provides 
a contribution oc AT 2 in the high temperature limit. We will see later that a term very 
similar to this can be extracted in the limit in which the distribution n k (0) is very large. 

Since the relevant scale is the quasiparticle mass Mr, we will choose to take M| > to 
describe an initial situation in which the O(N) symmetry is unbroken. 

It is convenient for numerical purposes to introduce the following dimensionless quanti- 
ties: 



M, 



T 



Mr t 



y, 2 (r) 



A 



R 



R 



2M 2 R 



8tt 



2 » 



M, 



R 



A 



2M 2 



(ir 2 ) R (t) - (tt%(0) 



(4.9) 



In terms of these dimensionless quantities the equations of motion (|4.1 , |4.2| ) with the initial 
conditions ((4.31) become 



^ + l + ^( r )_^(0) + ( 7 S(r) 



<p{r) = , 



(4.10) 



d 2 



+ g 2 + l + ^ 2 (r)-v? 2 (0)+^S(r) 




W q 2q 3 [ M 2 R 



q (r) = ; <p q (0) = 1 ; 0,(0) = -iW q , (4.11) 
lU , (4.12) 



where we have chosen the renormalization scale k = Mr for simplicity. 

In order to make our statements precise in the analysis of the spherically symmetric 
"tsunami", we will assume the initial particle distribution to be Gaussian and peaked at 
some value q and width £ so that 
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n q {(j) = — — exp 



Q-Qo 



(4.13) 



where N is the total number of particles in a correlation volume M^ 3 and J is a normalization 
factor . The case Nq >> 1 corresponds to the high density regime with many particles in 
the effective correlation volume. 



A. Preliminary considerations of the dynamics: 

Before engaging in a full numerical solution of the evolution equations, we can obtain 
a clear, qualitative understanding of the main features of the evolution by looking at the 
quantum mechanics of the wave functional. 

The dynamics is different for the different solutions for A q in Eq. fl3.22|) and can be 
understood with simple quantum mechanical arguments: 



1. Case I 



A 9 = 2K(0) + v /^(0) + ^(0)]: 
In this case 



A,(0) « 4 n g (0) » 1 for 



|g-go| 



for » i ; (4.14) 



and the covariance of the wave functional ( |3.1|) is given by 

AgiO) = Q q (0) = u q {l + A q ) » u q for ~ i 

~ Uq f or >> L ^ 15 ^ 

Now, for each wave vector q we have a Gaussian wave-function which is the ground 
state of a harmonic oscillator with frequency fi 9 (0), but whose evolution is determined by 
a Hamiltonian for a harmonic oscillator of frequency oo q (0) at very early times. For the 
modes q such that Q q » u q (0), the wave function is very narrow compared to the second 
derivative of the potential and there is a very small probability for sampling large amplitude 
field configurations (i.e. large rjf). It is a property of these squeezed states that whereas the 
wave-functional is narrowly localized in field space, it is a wide distribution in the canonical 
momentum (conjugate to the field) basis. This wave function will spread out under time 
evolution to obtain a width compatible with the frequency u q , i.e. the covariance A q {r) will 
diminish in time and the fluctuation 
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will increase in time. This will in turn cause the time dependent frequency in the Hamilto- 
nian, u) q (r) to increase with time since the frequency and the fluctuations are directly related 
via the self-consistency condition. The resulting dynamics is then expected to approach an 
oscillatory regime in which the width of the wavefunctional and the frequency of the har- 
monic oscillator are of the same order. Under these circumstances, there is the possibility 
that for a particular range of parameters (coupling, central momentum of the distribution 



and particle density) parametric amplification can occur fl39l , f44l that could result in particle 
production and redistribution of particles as will be discussed below within an early time 
analysis. We will see that this case corresponds to a "tsunami" configuration in a theory 
which is symmetric even in the absence of the medium. 

2. Case II 



A 9 = 2K(0)- v /n2(0)+n,(0)]: 
In this case 



A,(0) « -1 + J- for 



4^0 



for » i ; (4.17) 



and the covariance is 

*45o)<<^' ? ~ ?01 



u q for 19 ^ » 1. (4.18) 

Therefore, large amplitude field configurations with momenta in the narrow "tsunami" shell, 
now have a high probability of being realized. As before, the wave function for each k mode 
corresponds to the ground state of a harmonic oscillator of frequency 



n ff (o) 



u q 



An qo (Q) 



which evolves with a Hamiltonian for a harmonic oscillator with frequency ui q (0). In this 
case the wave function for q « q Q is spread out over field amplitudes much larger than 

and it is localized in the canonical momentum basis. 
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Under time evolution the wave function will tend to be squeezed i.e. it will be forced 
to diminish its width and to become localized inside the potential well. This implies that 
the covariance A q (r) will increase under time evolution, while the fluctuation ( |4.16| ) and the 
time dependent frequency oj q {r) will decrease, i.e. the potential 'flattens out'. 

In this case the quantity gSfV) (the renormalized quantum fluctuations) in the evolution 
equations ( ^4.11|) decreases and as will be seen below, under certain conditions, can become 
negative. 

There is thus a possibility of inducing spinodal instabilities in the quantum fluctuations. 
To see this, consider the case in which (p = and the effective mass squared A4 Rtt (t) = 



1 + (?£(t) in the equation ( 4.11]) becomes negative, i.e. when flfS(r) < —1. The modes for 
which q 2 < \M r (t)\ will see an inverted harmonic oscillator and they will begin to grow 
almost exponentially, resulting in copious particle production for these modes as can be seen 



from the expression for the particle number as a function of time ( |3 . 1 9 ) . 



This situation, in which the potential turns into a maximum at the origin dynamically, 
corresponds to symmetry breaking, since the minimum will be away from the origin. In 
this case the dynamics will result in a re-arrangement of the particle distribution: spinodal 
instabilities will arise, long-wavelength modes will begin to get populated at the expense of 
the initial non-equilibrium distribution. The spinodal instabilities will in turn result in an 
increase in the fluctuations that will tend to cancel the negative contribution to gll from 
the initial non-equilibrium distribution. Eventually a stationary regime should ensue in 
which the instabilities are turned-off and the distribution of particles will be peaked at low 
momenta. 

At this point we want to emphasize that the possibility for the onset of spinodal instabil- 
ities is purely dynamical. In contrast to previous studies of dynamics in spinodally unstable 
situations P^ ,p9| , P5|,p2| in which an initially symmetric state is evolved with a broken sym- 
metry Hamiltonian, in the present case the initial state and the effective Hamiltonian are 
symmetric and the instability is a consequence of the non-equilibrium dynamics. 

The above analysis of the dynamics, based on the quantum mechanical analogy will be 
shown to be accurate in the next section where we present the details of the numerical 
evolution. 

B. Early Time Analysis: 

A more quantitative understanding of these cases can be achieved by studying the early 
time behaviour of the solutions and setting ip = 0. 

1. Case I 

In this case with A q given by ( |4.14| )and focusing on the very early time during which 
backreaction effects can be ignored, the solution to the mode equations (Eq. ( [4. 1 1| ) ) with the 
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initial conditions given by Eq. 
^(0): 



is simply a superposition of plane waves with frequency 



T 



cos(o;g(0)r) — z(l + A q ) sin(a; g (0)r). 



(4.19) 



The renormalized quantum fluctuations which are dominated by the modes within the highly 
populated momentum shell are given by, 



0E(t) 



-g 



2 sin 2 K(0)r)[l-(l + A 9 ) 
Q q c (l + A ) 



(4.20) 



If the initial distribution of particles is sufficiently sharp, a qualitative understanding of the 
early time dynamics can be obtained by a saddle point analysis of the contribution from the 
region of large occupation number. 

In this limit <?S(t) is approximately given by, 



0E(t) 



4gN 



OJ, 



<1Q 



The mode equations now become 



i 2 

dr 



2gN 2gN 



., + q A + 1 + J " — - cos(2co>g r) 



'/" 



OJ, 



lit, 



<f> q (r) = 0. 



(4.21) 



(4.22) 



This is a Mathieu equation whose solutions are of the Floquet form |4"2 |. The first and 
broadest instability band is centered at the value of q given by 



q 



2 2gN 



OJ, 



The width of the unstable band depends on the parameter 

gN 



Q 



(4.23) 



(4.24) 



and can be read off in reference ||42|| . There is a rather small window of relevant parame- 
ters that could allow appreciable parametric amplification. Whether the backreaction effects 
allow the unstable band to remain under time evolution resulting in large particle produc- 
tion and redistribution of particles is a detailed dynamical question that will be studied 
numerically below. 



2. Case II 



The dynamics in this case can be unders tood by the he uristic arguments presented below. 
We work with the solution A q = 2[n q (0) — ^Jn^O) + n q (0)}, which in the limit iVo >> 1 yields 
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A q ^ -1 + for q w g (4.25) 

4n 90 

ps otherwise. (4.26) 
leading now to the following approximate form for the fluctuation at early times: 

gZ(r) ps -^sin 2 K r) (4.27) 



when the backreaction effects can be ignored. The first feature to note is that g and Nq 
appear together in such a way that the effective coupling is now gN and hence the physics 
is intrinsically non-perturbative when iVo ps 1/g. This situation is very similar to that in 
high temperature field theory wherein the relevant dimensionless quantity is T/m(T) (m(T) 
is the temperature corrected effective mass) and the effective coupling constant for long- 
wavelength physics is \T/m{T). In this situation the non-perturbative hard-thermal-loop 
resummation is required. 

Secondly the expression for gH(r) is always less than or equal to zero. Notice that unlike 
Case I, <?£(t) is negative [see Eq. ( f4.21|) 1. In particular g£(0) = and then <?£(r) becomes 
negative i.e. it begins to decrease. The fact that the fluctuations decrease was exactly what 
we had expected from the wave functional analysis presented in the previous section. 

Furthermore we see that when 4gN /u qo > 1 at very early times there will be an unstable 
band of wave-vectors. An estimate of the width of the band can be provided by averaging 
the time dependence of gH over one period of oscillation. This estimate yields the band of 
wave- vectors 



Q<q< W?±-l = qm (4.28) 

V u qo 

which will become spinodally unstable. The mode functions for these wavevectors will grow 
exponentially at early times and their contribution to the fluctuation g£(r) ( [4.9|) will grow 
- this is a back-reaction mechanism that will tend to shut-off the instabilities. 

This means that if we begin with a completely O(N) symmetric state i.e. <f{0) = y?(0) = 
and if we choose N large enough such that 1 + g£ ps 1 — 4giVo sin 2 (u qo r)/ij qo < 0, spinodal 
instabilities will be triggered and the symmetry will be spontaneouly broken. The condition 
for spinodal instabilities to appear is given by 

AgN 



> 1 (4.29) 



which determines the critical value of the particle number in a correlation volume in terms 
of the coupling and the peak momentum of the distribution. 

In the preceding sections we provided an intuitive understanding of the underlying mech- 
anism of symmetry breaking in terms of a quantum mechanical analogy. We now provide an 
alternative argument to clarify the physical mechanism for the dynamical symmetry break- 
ing. The argument begins with the expression for the 'dressed' mass in Eq. (|4.8| ) which we 
write in terms of dimensionless quantities as 
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q 2 dq 



mS = h4+ 9 mS/^ (0) 

2 2 , f q 2 dq 
™ R = m B + g J 



1 + A 



w g (0) ' 



(4.30) 
(4.31) 



The second term is dominated by the peak in the initial particle distribution. Using eqns. 
( Q3 , OBp and using a saddle point approximation assuming a sharp distribution, we obtain 
the relationship 



M 



R 



4gNo 

^qO 



m 



R- 



(4.32) 



Then choosing the effective mass M R > as we have done throughout, we see that when 
the condition for spinodal instabilities in Eq. ( |4.29|) is fulfilled then it must be that m 2 R < 0. 
Therefore the renormalized mass squared in the absence of the medium is negative and the 
medium effects, i.e. the non-equilibrium distribution of particles dresses this mass making 
the effective, medium 'dressed' mass squared positive. Thus in the absence of medium 
the potential was a (spontaneous) symmetry breaking potential. The initial distribution 
restores the symmetry at r = much in the same way as in finite temperature field theory 
at temperatures larger than the critical temperature. However the initial state is strongly 
out of equilibrium and its time evolution re-distributes the particles towards low momentum 
and the spinodal instabilities result from the squeezing of the quantum state as explained 
above. 

This situation must be contrasted to that in Case I above. The same argument, now 
applied to Case I leads to the result 



Ml 



1 + 



4:9 Np 

^qO 



m 



R- 



(4.33) 



We clearly see that with a positive effective mass, Case I corresponds to the situation in 
which the theory was symmetric even without the medium effects (i.e. m R > 0). 

Thus we obtain a physical picture of the different cases: in Case I the symmetry was un- 
broken without a medium and remains unbroken when the large density of particles is added. 
By contrast in Case II, the symmetry is spontaneously broken in the absence of a medium, 
the high density initial state restores the symmetry in a state out of equilibrium. Under 
time evolution the dynamics then redistributes the particles producing spinodal instabilities 
and breaking the symmetry. 

We reiterate that the second case represents a novel situation which is in a sense, contrary 
to what happens at high temperature where thermal fluctuations suppress the possibility of 
long-wavelength instabilities. 

The issue of symmetry breaking is a subtle one here. If we begin with symmetric initial 
conditions, ip(0) = 0(0) = 0, the wavefunctional will always be symmetric since the evolu- 
tion will maintain this symmetry. In order to test whether the symmetry is spontaneously 
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broken or not, one must provide an initial state that is slightly asymmetric, with a very 
small initial expectation value ip(0) 7^ 0, and follow the subsequent time evolution. If the 
expectation value oscillates around zero, then the symmetry is not spontaneously broken 
since the minimum of the 'dynamical effective potential' is at the origin in field space. If the 
expectation value begins rolling away from zero and reaches a stationary value away from 
zero then one can assert that there is a dynamical minimum away from the origin and the 
symmetry is spontaneously broken. Thus the test of symmetry breaking requires an initial 
condition with a small value of the order parameter. 



C. The Late Time Regime 



The asymptotic value of the order parameter can be obtained by analyzing the full 
dynamics of the theory and depends on the initial conditions. This reflects the fact that 
there is no static effective potential description of the physics. However, some information 
about the asymptotic state (when <p(oo) = <p(oo) = 0) can be obtained from the equation of 
motion ( 4.10|) by setting <p(oc) = which yields the sum rule |4^j45[ 



1 + <^ 2 (oo) + #£(oo) = 



(4.34) 



provided <f(oo) 7^ 0. This sum rule guarantees that the pions are the asymptotic massless 
Goldstone bosons since 



Mi 



2 , 

Tfl d H 



2 (oo) + (n 2 ) R (oo) 



(Eq. (|2.5| )) and the sum rule is a consequence of the Ward identities associated with the 
global O(N) symmetry. 

The non-linear evolution of the mode functions results in a redistribution of particles 
within the spinodally unstable band. The distribution becomes more peaked at low momen- 
tum and the effective potential flattens resulting in a non-perturbatively large distribution 
of Goldstone bosons at low momentum. 



V. NUMERICAL ANALYSIS 



1. Case I 



We have investigated the possibility of parametric amplification in this case in a wide 
region of parameters but always in the dense regime Nq » 1 and varying the center of 
the distribution. We find numerically that the backreaction effects shut off the parametric 
instabilities rather soon allowing only small particle production and redistribution of parti- 
cles. Typically the distribution develops peaks but remains qualitatively unchanged and the 
dynamics is purely oscillatory. 
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2. Case II 



The numerical analysis of the problem involves the solution of the coupled set of equations 
( |4.10|) , ( [4.11 ) and (|4.12| ) appended with initial conditions. We choose the zero mode initial 



conditions to be y?(0) = 10 3 ,0(O) = while the mode functions satisfy g (O) = 1; 4> q (0) 
-iuj q (0)(l + A q ). Here u q (0) = V<? 2 + 1 and 



A g = 2K(0) - y/n*(0)+n q (0)] . (5.1) 



We have tested the numerics with a momentum cutoff A = 25 in units of the renormalized 
mass Mr and found that after renormalization the numerical results are insensitive to the 
value of the cutoff provided it is chosen to be much larger than the largest wave- vector which 
becomes spinodally unstable. The initial particle distribution is chosen to be 

n q (0) = ^e~^\ (5.2) 

where the total initial number of particles is taken to be to be iVo = 2000, the coupling is 
fixed at g = 10~ 2 and the initial value of the order parameter is taken to be y?(0) = 10~ 3 . 
Results: 

Fig. flip shows g£(r) vs. r and the effective 'pion' mass squared M. 2 Rm {t). We see clearly 
that spinodal instabilities are produced and the quantitative features of the dynamics are 
in agreement with the estimates established for the early time dynamics given by Eq. (|4.27|) . 
We see that the pions become massless, asymptotically. The distribution function n q (r) 
multiplied by the coupling g is shown in Fig. @ at different times, clearly demonstrating 
how the distribution changes in time. As a consequence of the spinodal instabilities the long- 
wavelength modes grow exponentially and the ensuing particle production for these modes 
populates the band of unstable modes. In particular the amplitudes of the long wavelength 
modes that become spinodally unstable grow to be non-perturbatively large of order 1/g and 
dominate the dynamics completely. At earlier times the initial peak in the distribution at 
q pa 5 can still be seen, but at later times it is overwhelmed by the distribution at long 
wavelengths. Fig. (Q) shows a zoom-in of the distribution functions (gn q ) vs. q at r = 30, 80 
near q = and also near the peak of the initial distribution, around go ~ 5. We see a remnant 
of the original peak, slighly shifted to the right but much broader than the initial distribution 
and with about half the original amplitude. After r ~ 10 the distributions do not vary much 
in this region of momenta, but they do vary dramatically at low momenta. Fig. (|j) shows 
the total number of particles as a function of r. We see clearly that initially the total number 
of particles diminishes because the fluctuations decrease at early times. The long-wavelength 
modes begin to grow because of spinodal instabilities but their contributions are suppressed 
by phase space. Only when their amplitudes become non-perturbatively large, is the particle 
production at long-wavelengths an effective contribution to the total particle number. When 
this happens, there is an explosive burst of particle production following which the total 
number of particles remains fairly constant throughout the evolution. After the spinodal 
instabilities are shut-off, which for the values chosen for the numerical evolution correspond 
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to r m 2, the dynamics becomes non-linear. Whereas during the initial stages the dynamics 
is in the linear regime, after backreaction effects have shut-off the spinodal instabilities the 
further evolution of the distribution functions is a consequence of the non-linearities. 

Fig. (|5]) exhibits one of the clear signals of symmetry breaking. The order parameter 
begins very near the origin, but once the spinodal instabilities kick in, the origin becomes 
a maximum and the order parameter begins to roll away from it. Notice that the order 
parameter reaches a very large value, which is the dynamical turning point of the trajectory, 
before settling towards a non-zero value. We find that the value of the turning point and the 
final value of the order parameter depend on the initial conditions. To illustrate the non- 
perturbative growth of modes clearly, we have plotted the quantity g\4> q {r = 5)| 2 — g\4> q {j = 
0)| 2 in Fig. (H) which shows how the amplitude of the long wavelength modes becomes 
non-perturbatively large and of order 1/g. 

We have also carried out the numerical evolution with g = 10~ 2 , Nq = 4000 and g = 
10~ 3 , iVo = 40000 with the same value of go and found the same quantitative behavior, 
proving that the relevant combination is gN Q as revealed by the analytic estimates above. We 
have also confirmed that for gN << 1 there are no spinodal instabilities and the dynamics is 
purely oscillatory without a redistribution of the particles and with no appreciable particle 
production. When the peak of the initial distribution function is beyond the spinodally 
unstable band q > q m (see Eq. ( |4.28| )), the original distribution is depleted and broadened 
somewhat with irregularities and wiggles but remains qualitatively unchanged (see fig. |3|). 
However, when the peak of the initial distribution is within the spinodally unstable band there 
is a complete re-distribution of particles towards low momentum. The original distribution 
disappears under time evolution and after the spinodal time only the low momentum modes 
are populated. 



VI. SYMMETRY BREAKING, ENERGY, PRESSURE AND EQUATION OF 

STATE: 

A. Onset of Bose Condensation: 

We have seen both from the numerical evolution and from the argument based on the 
sum rule ( |4.34j ) which is a result of the Ward identities and Goldstone's theorem, that the 



effective mass term vanishes asymptotically. Therefore the asymptotic equation of motion 
for the mode functions is that of a massless free field. In particular the asymptotic solution 
for the q = mode is given by 

M r -> oo) = A + Br (6.1) 

where A and B are complex coefficients that can only be obtained from the full time evolu- 
tion. However because the Wronskian 

MrWo(r) - 4>o(t)Mt) = 2iW (6.2) 
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is constant in time, neither A nor B can vanish |43| . This situation must be contrasted with 
that for the q 7^ modes whose asymptotic behavior is of the form 



g (r -> 00) = a q e iqT + p q e~ iqT . (6.3) 

This causes the number of particles at zero momentum to grow asymptotically as r 2 whereas 
the number saturates for the q 7^ modes. The three dimensional phase space conspires to 
cancel the contribution from the q = mode to the total number of particles, energy and 
pressure, which, from the numerical evolution (see Fig. (4)) are seen to remain constant at 
long times. This situation is very similar to that in Bose-Einstein condensation where the 
excess number of particles at a fixed temperature goes into the condensate, while the total 
number of particles outside the condensate is fixed by the temperature and the chemical 
potential. The q = mode will become macroscopically occupied when r ~ \ V where V is 
the volume of the system (i.e. the number of particles in the zero momentum mode becomes 
of the order of the spatial volume). When this happens this mode must be isolated and 
studied separately from the q 7^ modes because its contribution to the momentum integral 
will be cancelled by the small phase space at small momentum. Again the situation is very 
similar to the case of the usual Bose-Einstein condensation. Notice that this argument is 
independent of a non-vanishing order parameter if and leads to the identification of the 
zero momentum mode as a Bose condensate that signals spontaneous symmetry breaking 
even when the order parameter remains zero. Since the effective mass is zero we identify the 
condensing quanta as pions and therefore this mechanism is a novel form of pion condensation 
in the absence of direct scattering. 

When scattering is included, beyond the leading order in the large N approximation, the 
formation of the Bose condensate will require a detailed understanding of the different time 
scales. The time scale for the collisionless process described above must be compared to the 
time scale for collisional processes that would tend to deplete the condensate. If spinodal 
instabilities causing non-perturbative particle production at low momentum occur on much 
shorter time scales than collisional redistribution then we would expect that there will be a 
non-perturbatively large population at low momenta that could be interpreted as a coherent 
condensate. 

The spinodal instabilities seen in this article are similar to those which lead to the for- 
mation of Disoriented Chiral Condensates [24.2^,25.3^]. However we emphasize that unlike 



most of the previously studied scenarios for DCC formation in which a 'quench' into the 
spinodal region was introduced ad-hoc, in the present situation the spinodal instabilities are 
of dynamical origin. We have studied a situation where the vacuum theory has symmetry 
breaking minima (with m 2 R < in Eq. fl4.32|) ) but the initial state is highly excited with 
the particle density larger than a critical value leading to a symmetry restored theory in the 
medium. However this initial state is strongly out of equilibrium and its dynamical evolution 
automatically induces spinodal instabilities. 
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B. Energy, Pressure and Equation of State: 



As mentioned in the introduction the goal of our study is to understand the dynamical 
evolution of strongly out of equilibrium states. In the usual investigations of the dynamics 
of the quark gluon plasma one uses a hydrodynamic description in which the energy density, 
pressure and all the thermodynamic variables depend only on proper time flU,!^]. The hy- 
drodynamic equations are then a consequence of the conservation laws which are appended 
with an equation of state to determine the evolution completely. The hydrodynamic regime 
corresponds to the case when the collisional mean free path is shorter than the wavelength 
of the hydrodynamic collective modes, and therefore the concept of local thermodynamic 
equilibrium is warranted. 

A valid question in the situation that we have studied in this aricle, is whether and when 
an equation of state is a meaningful concept. In the leading order in the large N expansion 
there are no collisional processes (these arise at 0(1/N)) and therefore the concept of a 
hydrodynamic regime in is not applicable in principle. Furthermore since the state considered 
is spatially homogeneous the pressure will depend on time rather than on proper time. Since 
the energy is conserved and the pressure evolves with time, an equation of state will have a 
meaning only when the evolution has reached the asymptotic regime. 

The energy density is given by 



E 1 12 , 1 2 ,2 , ,4 , 1 
= —Q> H TTLnO H W H 

NV 2 V 2 BV 8 v 4tt 2 



k 2 dk 
20, 



M 2 + w;(r)|^(r)| a -^(vrT. (6.4) 



The last term which arises in a consistent large N expansion, is extremely important in 
that after renormalization it provides a negative contribution which can interpreted as part 
of the effective potential [fPySQjl . Using the equation of motion for the order parameter and 
the mode functions it is straightforward to show that the energy is conserved and the last 
term is necessary to ensure energy conservation. Since the energy is conserved it can be 
renormalized by a subtraction at r = 0, and therefore is finite in terms of the renormalized 
quantities. 

The pressure is given by the following expression, 



p + E 
NV 



+ 



2vr 2 



k 2 dk 
2Q n 



(6.5) 



Unlike the energy density, the pressure is not a constant of the motion and needs proper 
subtractions to render it finite. The detailed expressions for both the renormalized energy 
and pressure can be found in references []3*9ifi5] , 4"B] . However, rather than computing the total 
energy density and pressure, we will study the contributions from the modes that are highly 
populated and whose amplitudes become non-perturbatively large (~ 1/g). Asymptotically, 
when the effective mass vanishes and the low momentum modes become highly populated 
with amplitudes of 0(1/ g) the renormalized energy density is given by (see ref. [0,fE|] for 
the explicit expression of the renormalized energy density) 
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E 



NV 4tt 2 



k 2 dk r 
~2Q n 



[r)\ 2 + ^ q {r)\Ur)\ 2 + 0{g) 



(6.6) 



where k m is the largest spinodally unstable wave vector at early times and O(g) represents 
terms that are perturbatively small. Using the asymptotic solutions for the mode functions 
given by Eq. (|6.3|) and neglecting the strongly oscillatory phases that average out at long 
times we obtain 



E M A R rim q ^dq 



NV 2ir 2 Jo 2W, 



+ m 2 +0(g). 



(6.7) 



Similarly, neglecting the contribution of modes with small amplitudes, we find that the 
renormalized pressure plus energy density is given by, 



p + E _ 4 M| /■*» q 2 dq 
NV 



3 2tt 2 Jo 2W 



so that in the asymptotic regime 



P 



E 
3~' 



(6.* 



(6.9) 



independent of the particle distribution which is non-thermal. This is one of the important 
results of this work. In Fig. (|7]) we show the trace of the energy momentum tensor E — 3P 
as a function of time, for the same value of parameters as for Figs. (1-5). Clearly, the trace 
vanishes asymptotically. 

During the early stages of the dynamics when spinodal instabilities arise and develop with 
profuse particle production, an equation of state cannot be defined. The dynamics cannot 
be described in terms of hydrodynamic evolution. Since the processes under consideration 
are collisionless, there is no local thermodynamic equilibrium and an equation of state is 
ill-defined. 



VII. CONCLUSIONS 

We have studied the evolution of an 0(iV)-symmetric quantum field theory, prepared in a 
strongly out-of-equilibrium initial state. The initial state was characterized by a particle dis- 
tribution localized in a thin spherical shell peaked about a non-zero momentum, a spherical 
"tsunami". The formulation of this scenario resulted from a simplification of the idealized 
'colliding-pancake' description of a heavy-ion collision. For a large density of particles in 
the initial state, the ensuing dynamics is non-perturbative and consequently we studied the 
O(N) theory in the leading order in the large N limit which is a systematic non-perturbative 
approximation scheme . When the tree-level theory has vacua that spontaneously break the 
symmetry and the number of particles within a correlation volume at t = is so high that the 
symmetry is restored initially , spinodal instabilities are then induced dynamically resulting 
in profuse particle production for low momenta. 
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This situation is to be contrasted with the usual studies of DCC's where the initial state 
is assumed to satisfy LTE (local thermodynamic equilibrium) and the spinodal instabilities 
are introduced either via an ad-hoc quench or via cooling due to hydrodynamic expansion 
which is also introduced phenomenologically. 

Backreaction of the long wavelength fluctuations eventually shuts off these instabilities 
and the nonlinearities redistribute particles towards low momenta. We thus find asymptoti- 
cally in time a novel form of pion condensation at low momentum, out of thermal equilibrium. 
Furthermore, a macroscopic condensate of the Bose-Einstein type will form at much longer 
times (mt ~ VV). 

When the spinodal instabilities shut off we find that the asymptotic 'quasiparticles' are 
massless pions with a non-thermal, non-perturbative distribution function peaked at low 
momentum but with an ultrarelativistic equation of state. 

We believe that these phenomena point out to very novel and non-perturbative mech- 
anisms for particle production and relaxation that are collisionless, strongly out of local 
thermodynamic equilibrium and cannot be described in the early stages via a coarse-grained 
hydrodynamic evolution. These are the result of strongly out of equilibrium initial states of 
high density that could potentially be of importance in the dynamics of heavy ion collisions 
at high luminosity accelerators. 

A more realistic treatment, modelling a collision will require an initial state which breaks 
the rotational invariance and selects out a beam-axis along which the colliding pions move 
in opposite directions. However, the analysis of such initial conditions is beyond the present 
numerical capabilities and will be deferred to a future work. 

An upshot of this study of high density, non-equilibrium particle distributions is the fol- 
lowing tantalizing theoretical question: can one extract a resummation scheme, or an effec- 
tive theory akin to the Hard Thermal Loop effective expansion for arbitrary non-equilibrium, 
non-thermal distributions such as the "tsunami" configuration for gauge theories'? Possible 
answers and consequences of such initial states for gauge theories will be discussed in a 



forthcoming article [47 
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FIG. 1. ^S(r) and M. 2 {t) vs. r respectively, with the initial distribution Eq.( |5.2D , 
N = 2000 ; g = 10~ 2 ; <p(0) = W' 3 ;q = 5 
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FIG. 3. gn q (r = 30,80) vs. q. with the same parameters as in Fig. (|T]), zoomed in the region 
near q ~ and the peak of the initial distribution go = 5. 
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FIG. 4. Total number of particles N(r) vs. r. with the same parameters as in Fig. (||) 
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FIG. 6. A g (r = 5) = g\4> q (r = 5)| 2 — g\<p g (r = 0)| 2 vs. q for the same parameters as in Fig. (||) 
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FIG. 7. E — 3P vs. r for the same parameters as in Fig. (|l|) 
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